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Abstract 

The energy-energy correlation function of the two-dimensional Ising model 
with weakly fluctuating random bonds is evaluated in the large scale limit. 
Two correlation lengths exist in contrast to one correlation length in the pure 
2D Ising model: one is finite whereas the other is divergent at the critical 
points. The corresponding exponent of the divergent correlation length is 
v e = 1/2 in contrast to the pure system where u e = 1. The calculation is 
based on a previously developed effective field theory for the energy density 
fluctuations. 

PACS numbers: 05.50.+q, 64.60.Cn, 64.60.Fr 



1 



Typeset using REVTfjX 



The two-dimensional Ising model |3J with random bonds is an interesting example for the 
competition of thermal fluctuations and quenched disorder. According to the Harris criterion 
neither of these two effects dominates over the other near the critical point. This picture 
is supported by a number of investigations for weak disorder, treated in perturbation theory 
with respect to the variance of the bond distribution, which came to the conclusion that the 
thermodynamic properties are not dramatically changed by quenched disorder Pfl. This 
is a consequence of the fact that the fixed point of the pure Ising model is stable against 
disorder (at least in one- loop order). The effect of disorder on the power law of the spin-spin 
correlation (SrSq) ~ Rr 1 ^ of the pure system is controversal. While some authors claim 
there is no effect on the exponent 1/4 others find a change for the correlation function 
to exp{— const. [log(logi?)] 2 } ||. On a naive basis one expects a reduction of the correlation 
near the critical point due to the additional fluctuations of the random bonds. Numerical 
simulations [|7||| do not indicate a significant change of the correlation function. 

Apart from the perturbative approaches it was found more recently in a non-perturbative 
treatment that disorder modifies the phase diagram of the two-dimensional Ising model 
||lJJ. In particular, disorder creates a new phase between the ferro- and the paramagnetic 



phase. However, the width of this phase is small and vanishes for any finite order in pertur- 
bation theory. For the same reason it is rather unlikely that this this phase can be observe 



in a numerical simulation. Nevertheless, there is a rigorous proof for its existence [IT]. This 
phase is characterized by a new order parameter. It was not possible in previous studies to 
understand its role in terms of the thermodynamic functions. However, it was related to the 



inverse correlation length of the average spin-spin correlation function by Braak |12| . In the 
following it will be shown that it also plays the role of an inverse correlation length of the 
average energy-energy correlation function. 

The Ising model describes a system of discrete spins S r = ±1 on a lattice. In this article 
a square lattice will be considered. The spins on nearest neighbor sites of the lattice are 
coupled with the energy 
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— — y ' Jr,r' S r S r i (1) 
r,r' 

where J r y is a random coupling variable. A statistical system of spins at the inverse tem- 
perature j3 is defined by the partition function 

Z = ]T exp(-/3W). (2) 

{5 r } 

Using this quantity the energy-energy correlation function (EECF) reads 

1 d d 

(£ r £ r >) = — ^ Jr,r"Jr',r"'-^ logZ. (3) 

P r n _/// OJ r r ii OJr'y" 



In order to calculate Z a number of mappings and approximations is useful. Firstly, the 
statistics of Ising spins can be mapped onto the statistics of closed polygons which excludes 
multiple occupation of links between lattice sites. This leads to Z — Z Zp where Z is 
a spin-independent part and Zp the partition function of the closed polygons. Secondly, 
Zp is equivalent with the partition function of non-interacting fermions with four colors. 
As a consequence, Zp is as fermion determinant. An approximation for weak disorder and 
large scales reads as a fermion determinant for complex fermions with two colors Zp m 
det(H + 5ya 3 ). The Fourier components of the lattice matrix H are 

H = ma 3 + y c (a 2 sin kx + a l sin k 2 ). (4) 

5y r represents the random fluctuations of the couplings with variance g and m = y — y c , 
y = tanh(/3J), y c = \[2 — 1. aj are the Pauli matrices. Introducing the Green's function 
G(e) = (H + 5ya 3 + ieao)^ 1 the EECF reads in this approximation 

(£ r £ r/ ) w Jim((y + 5y r )(y + 8y r ')Tr[G r yG r >, r \) (5) 

with a positive constant y. Naively, one would neglect in a weak disorder approximation the 
Sy contributions in front of the G"s. However, we will see in the following that these terms 
are important. This is the first step for the derivation of the EECF in the disordered Ising 
system. 
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As given in (|3[) the fluctuations of the random coupling J r y are conjugate to the energy 
density S r . Therefore, it is convenient to transform the distribution of the fluctuations 
Sy into a distribution of the energy fluctuations. The latter was performed in Ref. || by 
mapping the random bond fluctuations to a random matrix 

Sy r cr3 -> Qr- (6) 

Q r is a Hermitean 2x2 matrix with random matrix elements. It describes the energy density 
fluctuations. For instance, the average energy density is proportional to the average trace of 
Q. Without discussing the structure of this random matrix model (the details can be found 
in U) we turn directly to the saddle point approximation of the model. It is sufficient to 
consider the upper left 2x2 block matrix Q r = Qo + 8Q r where Q is given by a saddle 
point approximation || with 

i 1 
Qo = ~2 W3 ~ 2 msa °' 

rj and m s are solutions of the saddle point equations 

V = V9 1 (8) 
m 

m = m + m s = (9) 

with the integral 1 = 2 J[rh 2 + rj 2 + |/t| 2 ] _1 c/ 2 fc/ {2ir) 2 . The solution, which describes the new 
intermediate phase, is given by r] ^ as 

T] 2 = e-^/a - m 2 /4 (10) 

and m = m/2. rj vanishes at the critical points m c = ±2e _7r ^/9 anc i j s zero f or m 2 > m 2_ 

The EECF can be expressed by the Gaussian fluctuations around the saddle point solu- 
tion Qo as 

(£ r £ r ,) « jfj [yTr(Q r ) + Tr(Q 2 r )][yTr(Q r/ ) + Tr(Q 2 r/ )]exp[-(5Q, I5Q)]l[d5Q (11) 



with normalization M . The first order terms in Q are the contribution from yG whereas the 
quadratic terms Tr(Q 2 ) are a result of 5yG in (Q). Higher order terms in the exponential 
function are neglected in this approximation. Taking them into account (e.g. third order 
terms in Q) could lead to a renormalization of I. On the other hand, the saddle point 



approximation corresponds with the N — > oo limit of a more general model [ lOfl . Fluctuations 
appear as corrections with 1/iV -1 / 2 for each SQ. Therefore, expression ( |TTD should be a good 
approximation in terms of the 1/iV-expansion. 

The eigenvalues Ai,..., A3 of the 4x4 matrix I(k) are positive whereas A4 vanishes like 
rf if one approach the critical points m c : 

A 4 ~ const. \m + m c \ + bk 2 + o(k 3 ). (12) 

The mode of A4 is SQu — SQ22- Therefore, the terms with Tr{Q) = m s + 5Qn + 5Q22 are not 
critical in contrast to terms with Tr(Q 2 ). This fact explains why the 5y contributions are 
important in the expression of the EECF in ([H). The mode 5Qn +5Q22, which corresponds 

— 1/2 

to A3, has a characteristic finite correlation length £ 3 ~ const. A 3 . Furthermore, there is 
a divergent correlation length from the correlation of Tr(Q 2 ) which is proportional to \tj\~ 1 , 
the order parameter of the intermediate phase. It reads in terms of m £4 ~ £o| m + m c \~ l l 2 . 

Thus, the EECF is characterized by two different length scales. One is the finite cor- 
relation length £3, the other is the divergent correlation length £4 which is related to the 
order parameter of the new phase. In contrast to this, the pure system has only one 
correlation length which diverges at the critical point m = as £, pure ~ |m| _1 . 

In conclusion, quenched disorder reduces significantly the correlation length of the energy 
fluctuations in the two-dimensional Ising model (see table 1). This result indicates that 
disorder is not only marginal in this model. A similar reduction is expected for the spin- 
spin correlation function in opposition to the controversal results obtained previously in 



perturbation theory |pFpf. However, Braak [12| found in saddle point approximation that 



the correlation length of the average spin-spin correlation function in the intermediate phase 
is proportional to n~ 2 ; i.e., the critical exponent for this correlation length is v = 1 as in the 



pure Ising system. 
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TABLES 

TABLE I. Table of critical points and exponents for correlation lengths of the energy-energy 
correlation function (with the constant c = nyl « 0.54). 
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